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We show that the Landauer multi-terminal formula for the conductance of a nanoscale system is 
incomplete because it does not take into account many-body effects which cannot be treated as con- 
tributions to the single-particle transmission probabilities. We show that the physical origin of these 
effects is related to the viscous nature of the electron liquid, and develop a perturbative formalism, 
based on the time-dependent current-density-functional theory, for calculating the corrections to the 
resistance in terms of the "Kohn-Sham current distribution" and the exchange-correlation kernel. 
The difficulties that still remain in calculating the latter are critically discussed. 



I. INTRODUCTION 

The trend towards extreme miniaturization of elec- 
tronic devices provides strong motivation for theoretical 
studies aimed at characterizing and understanding the 
electrical transport properties of quantum-mechanical 
systems ii Here, by "quantum system" we mean a molecu- 
lar structure or cluster of atoms, or perhaps a microscopic 
region defined on the surface of a semiconductor. Either 
way, this system is connected to an external circuit which 
maintains current flow via electron sources. 

In the case of steady-state transport, this complicated 
non-equilibrium many-body problem is often times sim- 
plified by conceptually replacing the electron sources with 
ideal reservoirs, whose role is to define a local elec- 
tron distribution and a local electrochemical potential at 
which electrons are injected in, or extracted from the sys- 
temj2.i^ The reservoirs are conceptual constructs which 
allow us to map the transport problem onto an ideal sta- 
tionary scattering one, so that the time derivative of all 
local physical properties of the system and the current, 
is zeroj^ 

As a further simplification, one assumes that these 
reservoirs are adiabatically "connected" to leads in which 
non-interacting electrons are free to propagate before 
scattering at the lead-system interface.— The leads are 
only a convenient region of space where scattering states 
can be developed into an appropriate basis of the Hilbert 
space. This viewpoint to electrical conduction is known 
as Landauer approach. 

A schematic of this approach applied to a system con- 
nected to several leads is shown in Fig. [T] where the 
shaded region represents the system and the white re- 
gions are the leads, numbered 1 to N . The contacts be- 
tween the leads and the system can be very complicated, 
and should be considered part of the system. The proper 
lead, far from the contact, is a single-electron wave guide, 
which we can assume to have constant electrochemical 
potential /ii (i=l,..,N). At equilibrium all the leads are 
at the same electrochemical potential /z, and no current 
fiows in or out of the system. As we move slightly away 



FIG. 1: Schematic of a quantum system in a multi-terminal 
configuration. The central region is connected to leads of 
non-interacting electrons, in turn connected adiabatically to 
reservoirs of electrons. 



from equilibrium, the currents fiowing in the leads will 
be related to the electrochemical potentials by the linear 
relationships 



N 



(1) 



where the currents li are reckoned positive when they 
flow out of the system, and negative when they flow 
into the system. The coefficients Gij are the linear con- 
ductances of the system. In an ideal steady-state situ- 
ation [^i and li independent of time) the conservation 
of charge implies that the sum of all the currents is zero 
and therefore (for every terminal j) 
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Furthermore, the condition that the currents vanish when 
all the chemical potentials are equal implies that 



N 



Eg. 







(3) 



2 



It follows that the ofF-diagonal conductances Gy, with 
i ^ j are sufficient to completely characterize the linear 
steady-state response of the system. 

The Landauer multi-terminal formula offers an appeal- 
ing way to relate Gij to the quantum-mechanical prop- 
erties of the system. In this theory is proportional to 
the quantum-mechanical probability that a single elec- 
tron coming from lead i with energy E be transmitted 
into a different lead j at the same energy. In linear re- 
sponse and at zero temperature, this energy can be taken 
to be the Fermi energy Ep of the system. We call this 
coefficient Tij{Ep), and notice that, in general, it is a 
sum of all the partial probabilities of transmission from 
one of the momentum states of the incoming electron at 
energy Ep to one of the momentum states of the outgo- 
ing electron at the same energy (see, e.g., Ref. [Ij). Thus 
the Landauer formula reads 

G.,^'^T,,{Ep), [i^j) (4) 

It is important to note that the mathematical descrip- 
tion of this approach relies on scattering theory, namely 
on the transmission properties of single electrons in the 
leads that scatter at the leads-system interface. There- 
fore, for this description to be valid any interaction be- 
tween electrons can only he included at a mean-field level. 
Many-body interactions beyond mean field destroy the 
concept of single-particle transmission probability and, 
in fact, when taken into account also in the leads, they do 
not even allow for the derivation of a closed form for the 
total current. ^'^'^ All of the above issues are particularly 
relevant in nanoscale systems, where the current densi- 
ties at the junction can be substantially larger than in 
the bulk. A large current density implies a large number 
of scattering events per unit time and unit volume, thus 
making the description of transport phenomena in terms 
of non-interacting particle properties questionable ji 

It should thus not come as a surprise, and this is 
what we set to clearly show in this paper, that the Lan- 
dauer formula ^ represents an incomplete description 
of electrical transport in nanoscale systems. This point 
is particularly relevant these days, since there has been 
a surge of theoretical activities aimed at calculating the 
transport properties of these systems from "first princi- 
ples" . A popular way to tackle this problem is to extract 
the transmission function appearing in Eq. ([T]) from the 
one-electron Green's function,^ which is calculated from 
the self-consistent potential of the ground-state density- 
functional theory (DFT).^ In this manner, one hopes 
to include the most important effects of the electron- 
electron interaction without losing the simplicity of the 
single-particle theory. Indeed, one expects that interac- 
tions control the positioning of the single-particle energy 
levels of the system with respect to the Fermi level, and 
for this reason they have a large impact on the conduc- 
tance. 

However, from a more general theoretical standpoint 
things are not so simple. First of all, even if we assume 



that the physical approximations underlying the Lan- 
dauer formula ([1]) are a reasonable starting point to de- 
scribe electrical transport, the use of ground-state DFT 
in the present context is highly questionable, since one 
effectively uses a ground-state theory for an intrinsically 
non-equilibrium problem, even in linear response and in 
the dc limit .""i 

Indeed - and this leads us to the central message of 
this paper - it is precisely the non-equilibrium nature of 
the transport problem which renders Eq. ([4]) untenable. 
In a practical realization of a transport experiment, elec- 
trons are in a state of non-equilibrium and, therefore, 
their correlations are time dependent, even in the limit 
of zero frequency. These correlations give rise to scatter- 
ing processes that cannot be described by a mean-field 
theory and, under certain conditions, may infiuence sub- 
stantially their dynamics. 

Therefore, the Landauer formula, which has been de- 
rived within a single-particle framework, cannot be un- 
critically transferred to the many-body context, hoping 
that a proper inclusion of many-body effects in the single- 
particle energy levels will always suffice. In fact, in this 
paper we show that there are many-body corrections to 
the Landauer formula, which cannot be formulated in 
terms of single- particle transmission probabilities. 

In order to demonstrate this important point of prin- 
ciple we start from the rigorous formulation of the con- 
ductance in terms of the zero-frequency limit of the 
exact non-local conductivity tensor CTy (r, r';a;) of the 
interacting many-electron system (w is the frequency), 
in the linear response regime. We then resort to the 
time-dependent current-density functional theory (TD- 

CDFT)iiii^ii^ to show that the conductivity tensor a 
satisfies the integral equation 

(T = CTs - • • (5) 

where ct ^ is the resistivity tensor of a noninteracting sys- 
tem in the presence of a static potential Vg (also known 
as the Kohn-Sham potential) that reproduces the exact 

ground-state density, and P xc is a dynamical contribu- 
tion that will be defined precisely in Section [1111 

The linear response formulation of mesoscopic trans- 
port dates back to works by Fisher and Lee^ and 
Baranger and Stonei^ in the 1980s and has recently been 
combined with density functional theory by several sm- 
thors J-^iiSiiliSO, Because tTg is the conductivity of a non- 
interacting system, it is possible to analyze it microscop- 
ically by the method of Fisher and Lee^^ (later general- 
ized by Baranger and Stone^^), and thus show that this 
part of the conductivity alone leads to the Landauer for- 
mula ([!]) , with transmission probabilities computed from 
the Kohn-Sham potential Vg. This step is still within the 
assumptions of the Landauer approach, whereby the elec- 
tron sources are replaced by conceptual reservoirs whose 
role is to populate the single-particle states according to 
different Fermi functions, and these states can be de- 
veloped in terms of the single-particle states of the leads. 
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However, this is not the whole story, since there is also the 
contribution of the second term on the right hand side of 
Eq. ([5]). In other words, within the Landauer viewpoint 
to conduction, even if we knew the exact Kohn-Sham 
potential, including all the self-interaction and non-local 
corrections which are responsible for the correct align- 
ment of the one-electron energy levels, we would still be 
making an error in calculating the conductance from the 
Landauer formula 

Next, we examine the nature of the correction to the 
Landauer formula. We observe that Eq. ([5]) is alge- 
braically equivalent to the equation 

P = P., + ^xc (6) 

where P = (t is the exact non-local resistivity, Ps^cTg 

is the Kohn-Sham resistivity, and Pxc is the contribu- 
tion from many-body exchange and correlation. The re- 
sistivity controls the energy dissipation associated with 
a steady current distribution, and the presence of the 

xc correction P^c implies that there are mechanisms of 
dissipation that are not taken into account in the Lan- 
dauer approach of elastically scattering electrons, with 
relaxation and dephasing occurring only in the reservoirs. 
What Eq. ^ tells us is that electron-electron interactions 
make up for additional dissipation within the system, a 
dissipation that is physically a manifestation of electronic 
viscosity. As a matter of fact, the simplest approximation 

for Pxc, which is derived from the Vignale-Kohn approx- 
imation to TDCDFT,^^ is expressed precisely in terms of 
the viscosity of a homogeneous electron liquid: this ap- 
proximation shows that pxc is a positive kernel, always 
giving rise to a positive contribution to dissipation (i.e. 
an increase in resistance). 

The existence of viscosity contributions to the electri- 
cal resistance was first pointed out in Ref.[l^ where these 
contributions were called dynamical corrections because, 
as we have discussed above and will show below, they 
vanish in a strictly ground-state formulation of the the- 
ory. However, the relation of such contributions to the 
Landauer formula had remained somewhat unclear (see 
also Ref . [l7t ) . The present work shows conclusively that 
the Landauer formula ^ is incomplete, and the many- 
body corrections to it are precisely the "dynamical cor- 
rections" identified in Ref. 

The form of Eq. © suggests a simple perturbative ap- 
proach to the calculation of the resistances Rij (derived 
from the conductances Gy and defined more precisely 
below) , based on the minimal entropy-production princi- 
ple of linear-response theory. In brief, since the energy 
dissipation rate (proportional to the entropy production) 
computed from the single-particle (mean-field) theory is 
stationary with respect to a small variation of the Kohn- 
Sham current distribution js(r) (for given total currents 
in the leads) it follows that the additional dissipation due 
to the xc term is simply 

Wxc =js l^xc -Is , (7) 



to first order in Pxc- From this formula, and from the 
knowledge of the Kohn-Sham current distribution, we 
can straightforwardly extract the xc contribution to the 
resistances. The formula for the xc two-probe resistance 
of a quantum point contact or molecular junction which 
we presented in Ref. [H will be recovered as a special case 
of the general perturbative formulation. 

Finally, we consider some quantitative aspects of the 
theory. It must be said that a compelling comparison be- 
tween theory and experiments is still hampered in most 
cases by an imperfect characterization of the contact re- 
gion. Keeping this in mind, it is now accepted that the 
theoretical calculations of the conductance of molecular 
junctions, using the Landauer approach and ground-state 
DFT, overestimate the measured conductance by at least 
an order of magnitude.— Part of this discrepancy can cer- 
tainly be attributed to errors in determining the position 
of the energy levels of the system relative to the electro- 
chemical potential in the leads - errors which in turn are 
intimately connected to self-interaction corrections, dis- 
continuities in the xc potential as a function of particle 
number, and so on.^^ Even after correcting for these ef- 
fects, however, it seems that the computed conductance 
remains larger than the observed one, and it is here that 
our many-body corrections can play a decisive role. 

Our preliminary estimates of the size of the correction 
seem to indicate that the many-body viscous effects con- 
tribute only a small percentage to the total resistancef^^ 
For the case of two infinite jellium electrodes separated 
by a vacuum gap, the use of the viscosity as reported in 
Ref. [13 has shown an even smaller effect<^ But this does 
not mean that the issue is settled. 

First of all, it is important to note that these esti- 
mates have been based on an oversimplified description 
of the current density in nanoscale systems, by neglecting 
transverse variations of both the density and current den- 
sitjii^. For instance, as shown in Ref. transverse den- 
sity gradients increase the dynamical resistance. Quite 
generally, the transverse density and current density gra- 
dients and the spatial variation of the viscosity must all 
be taken into account when evaluating the viscous resis- 
tance. This is particularly relevant in nanoscale systems 
where non-linear (turbulent) effects have been recently 
predicted. Therefore, for a given nanoscale sys- 
tem, these dynamical effects need to be evaluated with 
the self-consistent microscopic density and current den- 
sity distribution. 

Aside from the above issues, there remains another and 
more fundamental source of uncertainty - namely, the 
value of the electronic viscosity which enters the dissipa- 
tive kernel pxc- In the concluding part of this paper we 
will argue that this value is still subject to a large uncer- 
tainty and we will outline the path along which better 
approximations might be obtained. 

This paper is organized as follows: Section ITTl reviews 
the general formulation for the conductance and the re- 
sistance of a nanoscale system in terms of nonlocal con- 
ductivity. In Section IIIII we present the time-dependent 
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current density functional approach to the calculation of 
the resistivity and demonstrate the existence of correc- 
tions to the Landauer formula In Section IIVI we 
develop the perturbative approach to the calculation of 
the many-body corrections to the resistance. In Sec- 
tion |V] we illustrate the working of the formalism in a 
simple one-dimensional model, re-deriving and extending 
the informal estimates of Ref. m. Finally, in Section IVTl 
we discuss the present difficulties in performing accurate 
calculations of the many-body corrections, and outline a 
path toward more accurate estimates. 

II. FORMULATION 

Our starting point is the linear response formula for 
the steady current density j in the presence of a steady 
electric field E: 

where a and fi denote cartesian indices and (Jap{Y, r') is 
the real part of the conductivity tensor. The integral runs 
over the whole volume of the system depicted in Fig. [TJ 
including the leads. The electric field, however, vanishes 
deep inside the leads. The above equation is satisfied for 
small electric fields of the form 

E(r) -Vr0(r) (9) 

where 0(r) is an electrostatic potential of arbitrary shape, 
except for the constraint of tending to constant values 

0(r) ^ /z, (10) 

deep into the i-th lead. 

Here, we assume that the electrostatic potential co- 
incides with the electrochemical potential deep into the 
leads. ''^ Because a steady current also satisfies the conti- 
nuity equation 

Vr-j(r)=0 (11) 

and because the current cannot be affected by a uniform 
shift of the electric potential in the whole space it follows 
that the conductivity tensor satisfies the conditionsi^ 

E5„5^a„0(r,r') =0 (12) 

a/3 

and 




where da is a short-hand notation for -^r- , and d'^ stands 
for Xj is the outwardly directed unit vector in lead 
J, and yj is a short-hand notation for the coordinates 



perpendicular to Xj, which are integrated over the cross 
section Cj of the j-th lead (see Fig. [T] for a schematic))^ 
The current li in the i-th lead is given by 




(14) 

Following Baranger and Stone^'^^^ we make use of Eqs. (fTO)) 
and p2|) . and an integration by parts to find the intuitive 
result (cf. Eq. U)) 

N 

U = Y,G,,li, (15) 

where 




Note that up to this point we have made no approxi- 
mation on the microscopic physical mechanisms that con- 
tribute to the conductance (jl6p apart from those embod- 
ied in the viewpoint represented in Fig. [TJ Therefore, 
within this viewpoint, the conductance (jl6p contains, in 
principle, all many-body interactions, even beyond mean 
field. 

The next step is then to express the conductivity ten- 
sor in terms of a microscopic current-current response 
function. To this end we introduce the proper current- 
current response function, which yields the electric cur- 
rent response to the fully screened vector potential in the 
following manner 

j„(r,w) = -e^^ / dr'xa/3(r,r';w) [Ap{y')^AhA^')] : 

(17) 

where A(r) is the external vector potential and A// (r) is 
the vector potential additionally created by the screening 
charge;^ The factor (e being the absolute value of the 
electron charge) is introduced to be consistent with the 
definitions used in other publicationsiii^i^. Then the 
conductivity is 

CTa/3(r, r ) = -e^ lim ^ . (18) 

The proper current-current response function is best 
expressed in terms of an infinite series diagrams with two 
current vertices, such as the diagrams shown in Fig. [2] 
where the solid lines represent free particle propagation 
and dotted lines Coulomb interactions. Notice that this 
series does not contain any diagrams that can be divided 
into two parts by cutting a single Coulomb interaction 
line. It is the exclusion of these diagrams that makes our 
response function "proper" , as opposed to "full" M- 
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FIG. 2: Diagrams for the proper current-current response 
function. The solid dots represent (particle) current vertices. 
The solid lines represent the free particle propagator and the 
dotted lines Coulomb scattering processes. 



A. Mean-field approximation 

In the special case of a non-interacting system, or a sys- 
tem interacting at a mean-field level, only the first term 
of the series survives and we get, following the standard 
rules^ 

nm 

(19) 

where n and m denote exact single-particle eigenstates 
with energies e„ and wave functions 'i/'n(r), f{^n) is the 
Fermi distribution at the common chemical potential 
/X (before applying the bias) and temperature T, and 
(i") is the matrix element of the a-component of the 
(particle) current operator between states m and n: 

m 

(20) 

Eq. leads upon substitution in Eq. to the 
standard Landauer formula ([4]).— The calculation is 
quite subtle, hinging on the possibility of choosing a com- 
plete set of exact eigenstates in the form of scattering 
states, i.e. states of energy e which describe a single par- 
ticle "entering" the system in the transverse channel a of 
the i-the lead, and scattered with probability amplitude 
Uajb into any transverse channel b of the j-th lead. 

Within this mathematical assumption, the transmis- 
sion coefficient that appears in Eq. Q is found to be 
given by 

Tij = ^ \tia,jbf ■ (21) 
a,b 

We refer to the original papersi^ii^ for the details of this 
derivation. What is important for our purposes is that 
the conventional Landauer multi-terminal formula ^ 
emerges from an approximation to the exact formula (jl6[) 
- an approximation in which only the first term in the 
infinite series of diagrams for the proper current-current 
response function is retained. 

The question now arises how to go beyond this simplest 
approximation to include electron-electron interaction ef- 
fects. In the next section we describe an approach based 
on time-dependent current-density functional theory. 



III. TIME-DEPENDENT CURRENT-DENSITY 
FUNCTIONAL THEORY 

As discussed in the Introduction, a popular approach 
to the inclusion of many-body effects in nanoscopic trans- 
port is to use the Landauer formula (|4]), but calculate 
the transmission probabilities by solving the one-particle 
scattering problem in a static effective potential that in- 
cludes many-body effects. How is such a potential to be 
constructed? 

The ground-state density functional theory (DFT) of 
Hohenberg, Kohn and Sham offers a practical answer 
According to this theory it is possible to find, in princi- 
ple, an exchange-correlation potential which, in combi- 
nation with the Hartree potential and the external po- 
tential, produces the correct ground-state density of the 
many-body system. Furthermore, this potential (known 
as the Kohn-Sham potential) is uniquely determined by 
the density. Thus, it is very tempting to make use of the 
Kohn-Sham potential to calculate the transmission prob- 
abilities and hope that all many-body effects pertaining 
to the transport problem be included. Unfortunately this 
approach lacks any rigorous theoretical foundation. In 
practice, it amounts to "dressing up" the free particle 
lines in the first diagram of Fig. [2l while still discarding 
all the other diagrams. Therefore, it must be interpreted 
as nothing more than a single-particle mean-field approx- 
imation even if we knew the exact ground-state xc func- 
tional, and as such there is no physical reason why this 
should be even approximately correct. 

The time-dependent current density functional theory 
offers a more solid basis to attack the problem. Taking 
for granted the ordinary DFT description of the ground- 
state, the TDCDFT attempts to describe the current 
response of the many-body system as the response of 
a non-interacting reference system to an effective time- 
dependent vector potential. The non-interacting refer- 
ence system is usually taken to be the "Kohn-Sham sys- 
tem" , i.e. the non-interacting system that is used in or- 
dinary DFT to reproduce the ground-state density of the 
many-body system. Thus, in the TDCDFT approach 
the current response to a time-periodic vector potential 
A(r,t) = A(r, aj)e~*"* + c.c, is written as 



ja{r,uj) = - 




+ AH,fj{r',u;) + A,cA^',io)} (22) 

where Xs,c(/3(i"j i"'; is the current response function 
of the Kohn-Sham system. Ah is the Hartree vector 
potential, and A^c is the exchange-correlation vector 
potential)^ The essential point is that the exchange- 
correlation potential is a unique functional of the current 
density, and in the linear approximation can be repre- 
sented as 

- e^A^cai'^, ^) ^ dr' ^ fxc.apir, r'; w)j^(r', u) , 

/3 

(23) 
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where the "exchange-correlation kernel" fxc.ai3{Y,r';uj) 
is determined by the ground-state density. Later we will 
also need the "exchange-correlation electric field" , which 
is defined as 



(24) 



The above formula should be compared with the 
exact linear response formula ^T7\ . Combining Eqs. (P^ . 
((23)) and p?)) we find the well-known relation between x, 
Xs, and fxc, namely 

[X^^]apir,r',uj) = [xs]~^{r,r', uj) - /:rc,a/3(r, r', w) (25) 

where [x~^]ai3iY,r' ,uj) is the matrix inverse of 
XQ/3(r, r', (jj), which is regarded as a matrix with 
indices a, r and /3,r' . 

Eq. ([25]) gives us a handle on the inverse of the conduc- 
tivity, i.e., the resistivity tensor. To make the connection, 
observe that the complex conductivity tensor d-{r,r',uj) 
is given by 



o-Q/3(r,r',Lj) = 



,Xap{r,r';uj) 



(26) 



whose real part reduces to (j{r, r') in the limit w — > 0. 
Accordingly, the complex resistivity tensor is given by 

Pa/3(r,r',cj) = [a^'^]a0{r,r',uj) = -^[x"i]a^(r, r', w). 



Then making use of Eq. ([25]) we find 



(27) 



Pa/3(r, r', uj) = Ps,afj{r, r', uj) + -^fxc,ai3ir, r', uj) (28) 

where the "Kohn-Sham resistivity", ps, has the same re- 
lation to Xs as the full interacting resistivity to x- Fi" 
nally, taking the real part of both sides and going to the 
zero-frequency limit we find 



Paf3{r, r') = Ps,Q/3(r, r') + Pxc,ai3{r, r') 



where 



Pxc,ai3ir,r') = - lim -ir'^mfxc,ap{r,r' ,uj) 



e 



(29) 



(30) 



This is the main result of our paper which we have 
anticipated in Eq. We clearly see that the resistiv- 
ity of the Kohn-Sham system - a non-interacting sys- 
tem in which many-body effects enter only implicitly 
through the static exchange-correlation potential - is not 
the whole story. This means, in particular, that it is 
not possible to give an exact representation of the con- 
ductance in terms of single-particle transmission proba- 
bilities. We have demonstrated this point for the linear- 
response regime, namely in the limit of zero external bias. 
However, this result must be valid also out of linear re- 
sponse, even though in this case the extent of the dynam- 
ical corrections, which in linear response are embodied in 
fxcap, is not so easily determined. 



We now shift our attention to the estimate of the dy- 
namical exchange-correlation contribution pxc which is 
controlled entirely and explicitly by many-body effects, 
i.e., time-dependent correlations in the effective poten- 
tial of TDCDFT. The many-body kernel fxc^ which ap- 
pears in Eq. (j29p . is not known exactly for any system, 
but a local approximation to it is available and has been 
used in the recent literature with varying degrees of suc- 
cess j^^^i^i^i^ In a local approximation the key quan- 
tity 

i?;rc,a(r,w) = " ^ X! / '^'''/^^ca/S (r , r' , w)j/3 (r' , w) , 

(31) 

which has the physical significance of "exchange- 
correlation electric field" (see Eq. ([24]) ). is taken to be 
a function of the local value of j(r) and its first and 
second spatial derivatives. The simplest approximation 
in this class is the so-called adiabatic local density ap- 
proximation (ALDA), which provides an instantaneous 
connection between E2.c(r) and j(r). In this approxima- 
tion, however, fxc is purely real, resulting in an exchange- 
correlation electric field that is always 90° out of phase 
with the current density. Therefore, such a field cannot 
contribute to the d.c. resistivity, consistent with the fact 
that 3to/xc vanishes in this approximation. 

So in order to obtain exchange-correlation corrections 
to the resistivity we must go beyond the adiabatic ap- 
proximation. This can be done with the help of the 
VK local density approximation,— which calls into play 
the viscosity of the electron gas. In this approximation 
the exchange-correlation field has a dissipative compo- 
nent which is ISO*' out of phase, i.e., opposite to the 
current. In the zero- frequency limit this component of 
the exchange-correlation field has the forniH 

Exc,a{r) = - ^ J dr'pxc,a0{r,r')jp{r') 



e^n(r) 



/3 



Ja(r) 
n(r) 



da 



n(r) 



3^^ • [mt] ' 



(32) 



where j(r) is the electric current density, and 77 (r) is the 
d.c. shear viscosity of a homogeneous electron gas of den- 
sity n(r).?? Eqs. ^9^, and constitute a complete 
(albeit approximate) formulation of the microscopic re- 
sistivity tensor within TDCDFT. The Kohn-Sham resis- 
tivity itself is accessible from the ordinary static DFT. 
Pxc is best described through the effective electric field 
it produces - an electric field E^jc directed against the 
current, which therefore does negative work on the cur- 
rent. In the next section we show how our expression 
for Exc can be directly applied to the calculation of the 
macroscopic conductance. 
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IV. PERTURBATIVE CALCULATION OF THE 
CONDUCTANCE 

Let us return to the system shown in Fig. [1] and let us 
assume that the electrochemical potentials /^i are peri- 
odically modulated in time with a (very) small angular 
frequency u. By "small" we mean a frequency much 
smaller than any other internal frequency of the system. 
The lead currents induced by the modulation are then 
also periodic and given by Eq. ([1]). Since the reservoirs 
are the only part of the system on which we have direct 
control it is evident that the work done on the system 
per unit time is 



N 



(33) 



where the angular brackets denote a time-average over a 
period of oscillation. This is also the energy that must 
be internally dissipated if the system is to remain in the 
steady state. 

In order to express W in terms of the lead currents we 
must invert the linear relation ([T]) between the lead cur- 
rents and the electrochemical potentials. Strictly speak- 
ing, this relation is not invertible, because a rigid shift 
of all the electrochemical potentials has no effect on the 
current. But the problem is easily solved by permanently 
grounding one of the reservoirs, say the one with i = 1 
so that /ii = at all times. Then the linear relation be- 
tween the remaining A'^ — 1 currents, /2, --.jIn and the 
corresponding electrochemical potentials /i2, ■■■^fJ-N "is in- 
vertible, and the current in the grounded lead is simply 
given by Ii = —I2 — ■■■ — In- 

Then, we see that the dissipated power can be repre- 
sented as 



N 



(34) 



where the {N — 1) x {N — 1) matrix Rij is the inverse 
of the matrix Gy stripped of the first row and the first 
columnj2I The macroscopic expression for W in terms 
of Rij is now equated to the usual microscopic expres- 
sion in terms of the resistivity, resulting in the following 
equation: 

where the integrals run over the volume of the system, 
including the leads, and we have dropped the time av- 
erage by going to the zero-frequency limit. Finally, by 
making use of Eq. (|29p we arrive at 

(36) 



This equation is formally exact if we know the exact 
many-body kernel fxc- Let us compare it with the for- 
mula we would obtain from the conventional single par- 
ticle theory, i.e. from the Landauer formula ([3]), for the 
same lead currents. 

As discussed in Section the "single-particle" theory 
assumes that all many-body effects can be included in the 
Kohn-Sham potential of the ground-state. Apart from 
this, the system is non-interacting. The current density 
distribution js(r) of this fictitious Kohn-Sham system is 
in general different from the true current density distri- 
bution j(r) of the many-body system, even though the 
macroscopic lead currents, i.e. the fluxes of j and js into 
the leads, are imposed to be the same. Therefore we 
write 



N 

E 



Ii Rs , ij Ij 



dv / dr'j,(r). Ps (r,r')-j.(r'), (37) 



where Rs,ij are the macroscopic resistances of the Kohn- 
Sham system, obtained from the standard single-particle 

theory. Notice that the dynamical term P xc is absent in 
the single-particle theory. 

The comparison between Eqs. and ([57)1 is com- 
plicated, in general, by the difference between j and js. 
A simple comparison becomes possible in the perturba- 
tivc limit, i.e. under the assumption that the dynamical 

many-body correction embodied in P xc is small. To this 
end we observe two facts: (i) The difference between j 
and js is of first order in the xc correction, and (ii) the 
right hand side of Eq. ([37]) is stationary under a small 
variation of the current distribution, such as the differ- 
ence between j and jg. The physical reason for this is 
that in any linear system with external leads (such as the 
Kohn-Sham system we are considering here) the power 
dissipated is stationary with respect to an infinitesimal 
variation of the current distribution at constant lead cur- 
rents. This implies that a first-order variation in the cur- 
rent distribution (about the steady distribution js in this 
case) produces a second-order variation in the dissipated 
power. Taking this into account we see that we can safely 
replace js by j in Eq. ([57)1 and subtracting from Eq. 
we arrive at the main result of this section: 

h^R^jIj = J drj dr'js(r)- (r,r') •js(r') 

i,j=2 

where 



(38) 



AR, 



(39) 



is the dynamical correction to the resistance. This equa- 
tion expresses the many-body correction to the macro- 
scopic resistance in terms of two things that are approx- 
imately known and/or calculable, namely, the resistivity 

exchange-correlation kernel Pxc: defined by Eq. ([5^ . and 
the Kohn-Sham current distribution of the device, js(r). 
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n(x) 



FIG. 3: A model two-terminal device in which the density 
changes only in one dimension. 1 and 2 are the lead regions. 



associated with the macroscopic lead currents 7^. The 
latter can be calculated, in principle, from the response 
of a non-interacting system to a screened electric field. In 
practice, one can calculate the correction to Rij by con- 
sidering a special situation in which only li and Ij arc 



different from zero. Then the left hand side of Eq. ((38l 
gives us exactly the desired correction to Rij . 



V. AN EXAMPLE 

Let us consider a simple model application of the gen- 
eral formalism. Our system is a potential barrier con- 
nected by two identical homogeneous leads (labeled 1 and 
2 in Fig. [3]) to two reservoirs, a "source" and a "drain", 
aligned along the x axis. The system is perfectly homoge- 
neous in the transverse directions y and z. The density 
changes only in the x direction. The source (terminal 
1) is grounded, so we only need to determine the two- 
terminal resistance R22 or, equivalently, the conductance 
G22 = R22 ■ Let / be the current through the device and 
i{x) = 1 1 A the current density, where A is the transverse 
cross section of the device. 

Notice that in this case there is no question of j(a;) be- 
ing different from (a;) since, by continuity, they are both 
uniform and equal to I/V. In the absence of electron- 
electron interactions the conductance of this system is 
simply given by Eq. ^ - the total transmission proba- 
bility across the potential barrier being the sum of the 
transmission probabilities of all the occupied transverse 
modes. Including the electron-electron interaction has 
two effects. The first is fairly trivial, namely, the effective 
potential in the ground-state is modified by screening and 
exchange-correlation effects, and the transmission prob- 
abilities must be recalculated for this effective potential. 
Up to this point the single-particle formula ([4]) remains 
in force. 

The second effect is the dynamical exchange- 
correlation correction - an effect that cannot be forced 
into the mold of the static mean-field theory. Making use 
of Eq. dSSl) with N = 2 and i{x) = I/Awe obtain 



Ai?2 



1 



dv I dv' p.j,c{x,x') 



(40) 



where pxc{x, x') denotes the xx component of the tensor 
P xc {x,x'), and we emphasize the fact that, due to our 



assumptions, it depends only on x and x' . Now observe 
that, according to Eq. p2p the action of Pxc on a uniform 
current density is specified by 



dv' pxc{x,x') = - 



3e'^n{x) 



dx 



ri{x)dx 



1 



n(x) 



(41) 



where ri{x) is the shear viscosity of the homogenous elec- 
tron gas evaluated at the ground-state density n{x). Sub- 
stituting this in Eq. (|40p and doing an integration by 
parts we arrive at 



AR22 = 



3eM 



dxrj{x] 



HxW ' 



(42) 



where n'{x) — dxn{x). Notice that this has the correct 
dimensions (Ohm) because n is a three-dimensional den- 
sity and T] has the dimensions of h times a density. Fig. [3] 
shows the electronic density in the leads and in the de- 
vice. Clearly the dynamical correction comes entirely 
from the non-homogeneous regions near the edges of the 
barrier (the contacts). 

Eq. ([^2]) was first obtained in Ref. [H, in a more intu- 
itive manner. The advantage of the present formulation 
is that it allows easy extension to more complicated sit- 
uations. For example, we can include the dependence of 
the density on the transverse coordinates y and z, while 
neglecting variations of the transverse components of the 
current. In this case, we still have j(r) — constant, but 
now the gradient of the density has both longitudinal and 
transverse components. As a result we get 



AR22 = 



2^2 



J dr r]{r) I ^ 



|Vp(r) 
[n(r)]4 



|V±n(r)| 
[n(r)]4 



(43) 

where V|| is the x-component of the gradient and Vj^ 
is the gradient in the y — z plane. This result (for con- 
stant viscosity) was first reported in Ref. [2^ following an 
intuitive procedure, still based on the calculation of the 
power dissipated in the circuit. 



VI. DISCUSSION AND CRITIQUE 

In this paper we have shown that the single-particle 
mean-field framework of Landauer is inadequate in prin- 
ciple to describe the transport problem in nanoscale sys- 
tems. That is to say a calculation of conductance from 
Eq. (HI) would not provide the exact current even if one 
could determine the transmission probabilities with the 
utmost precision. Dynamical, many-body effects enter 
the picture due to the intrinsic non-equilibrium nature 
of conduction.^ These effects cannot be captured by a 
static formulation. 

The next question is: what is the actual size of these 
dynamical corrections? One of the main results of this 
paper, Eq. (j38p . opens the way to a fully microscopic first- 
principles calculation of nanoscopic resistances and con- 
ductances (in the linear regime) within the framework of 
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FIG. 4: Qualitative behavior of the viscosity of a homoge- 
neous electron e function of frequency. For frequencies 
smaller than 1/r, with r the quasiparticle relaxation time, rj 
approaches the d.c. limit rjo- r/o vs. density is shown in inset 
(a) for T — 300 K. The calculation is done with the Abrikosov- 
Khalatnikov formula^^. For frequencies larger than 1/r, but 
still much smaller than the Fermi frequency Ep/h, rj tends 
to a different constant 7700. The behavior of 7700 vs density is 
shown in inset (b) in various approximations at zero temper- 
ature. CV is from Ref.— , QV is from Ref. [3^ . The atomic 
unit of viscosity is h/a% =~ 7 x 10^'^ Poise (ob is the Bohr 
radius) and na% — j;^. 



frequency. In other words, the temperature (and, with it, 
the quasiparticle scattering rate) went to zero before the 
frequency. When calculated in this manner, the viscosity 
turns out to be very small indeed: its value is in the 

range of 10~^/r^ Poise, where ( i^F^^j — 1 is 

the average inter-electron distance in units of the Bohr 
radius as ■ (For comparison, water at room temperature 
has a viscosity of about 10~^ Poise). 

On the other hand, it is well known from the theory 
of homogeneous Fermi liquids22,, that the behavior of the 
viscosity is quite different if the zero-frequency limit is 
approached at finite temperature. Namely, in this case 
the viscosity turns out to be proportional to the mean 
free path of the quasiparticles, which grows as l/T^ in 
the low temperature limit. The divergence of the zero- 
frequency viscosity for T ^ reflects the fact that long- 
lived quasiparticles can transport momentum arbitrarily 
far away from the source of the stress. This is also the 
reason why the viscosity of an ideal classical gas is inde- 
pendent of density;^ since an increase or a decrease in 
the frequency of molecular collisions is exactly compen- 
sated by an opposite variation in the molecular mean free 
path. 

Panel (b) of Fig. [3] shows the behavior of the zero- 
frequency viscosity at room temperature for an electron 
gas, estimated from Eq. (7.22) of Abrikosov and Khalat- 
nikov (AK)f^ with due allowance made for the different 
form of the interaction potential (the AK work was for 
^He). The Eq. (7.22) of AK can be rewritten as 



the local approximation to time -dep endent current den- 
sity functional theory. In Ref. we tried to address 
this question for the simple quasi-one-dimensional model 
discussed in the previous section and found that the vis- 
cosity correction to the resistance was only a small frac- 
tion of the total. For the case of two infinite jellium 
electrodes separated by a vacuum gap, a more accurate 
calculation based on the homogeneous electron gas vis- 
cosity reported in Ref. [13 (see Fig. [4|) has shown an even 
smaller effect»^ But, as pointed out in the introduction, 
this does not mean that the issue is settled. 

Looking back at Eq. ([55]) we see that an accurate eval- 
uation of the viscosity correction has two ingredients: (1) 
the Kohn-Sham current density distribution and (2) the 
viscosity of the homogeneous electron gas. As for the cur- 
rent density, it is important to note that all the estimates 
so far have been based on an oversimplified model in 
which the current density was assumed to be uniform in 
spac o^^'^^1^^ . In general, the spatial variation of the cur- 
rent density cannot be neglected; especially in nanoscale 
systems where large transverse variations of the current 
density are commonj^i^i2£i2£ 

Another and more fundamental source of uncertainty 
is in the value of the electronic viscosity which enters 
the dissipative kernel pxc- The viscosity we have used so 
far, which is plotted in inset (b) of Fig. [H was obtained 
from a zero-temperature calculation in the limit of zero 



770 = hn 



15tt \kBT 
'w(C 



[kpaB) 



cos ■ 



(1 — cos( 



)2sin2, 



(44) 



which shows explicitly the physical dimensions of the vis- 

cosity (hn). Here kp is the Fermi wave vector, Ep = 2m 
is the Fermi energy, is the Bohr radius and w(6, (j)) is 
the square of the matrix element of the electron-electron 
interaction potential (expressed in units of 47re^/fc|.) be- 
tween the initial and the final state of a collision process 
with incoming momenta Pi,P2 and outgoing momenta 
p']^,P2, where 9 is the angle between the incoming mo- 
menta and (j) is the angle between the planes formed by 
(Pi;P2) and (pi,P2). All the momenta are close to the 
Fermi surface, and the symbol av denotes the average 
over 6 and cj). The simplest approximation for w is the 
Thomas-Fermi approximation, in which we have 



1 



2(1 — cos 0) -I- 4ars/7ry 



(45) 



with a = (4/97r)^/^ ~ 0.521. With this approximation 
the average over 9 can be done analytically and Eq 
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can be rewritten as 



hn 




1.579 X 10^ 
T 



\ 2 TT 



(46) 



The result of the evaluation of this expression is shown 
in inset (a) of Fig. |4] for T — 300 K (room temperature). 
Notice that the presence of the factor 1/rf causes the 
viscosity to increase sharply with increasing density, in 
contrast with what we observe in inset (b) of Fig. U) It 
is evident that the d.c. viscosity is orders of magnitude 
larger than the finite frequency viscosity plotted in inset 
(b) of Fig. H 

What does this imply for our analysis of the conduc- 
tance in nanoscopic and mesoscopic devices? Obviously, 
these systems do not host a uniform electron liquid, and 
in particular they do not support long-lived quasipar- 
ticles that can transport momentum to infinity. This 
means that the large finite-temperature results of the 
uniform electron liquid are almost certainly not relevant 
for nanoscopic devices: the mean free path of quasipar- 
ticles is naturally limited by the geometric size of the 
device. "'^'^^ However, the huge difference between the nu- 
merical values of the viscosities in the insets of Fig. ^ 
suggests the possibility of a mesoscopic "middle-ground" 
which under certain conditions may be much larger than 
the zero-temperature viscosity. This, however, is unlikely 



to be "universal" , rather it must be related to the specific 
microscopic geometry of the system. A related difficulty 
is that, in general, the rate of dissipation in an interacting 
electron system depends strongly on the excitation spec- 
trum of the system. Modeling dissipation through the 
viscosity of a homogeneous electron g implied by 

our local-density approximation, may lead to a severely 
distorted description of the dissipative process. A truly 
universal description of dissipation (if possible at all) is 
still out of sight. 

A central issue emerges from the above discussion, 
namely the need for an accurate, testable, and reliable 
dissipative functional for time-dependent current density 
functional theory. The local-density approximation is 
only a first step. However, once a better functional is 
proposed, our formalism provides a simple and elegant 
way to test its predictions for the resistance of nanoscale 
systems. 

Finally, we stress that we have focused our attention 
to the linear-response regime. It would be interesting 
and important (although not trivial) to extend the re- 
sults presented in this paper to the non-linear case. Such 
an extension would allow analysis of the many-body cor- 
rections to the current- voltage characteristics ~ and cor- 
responding dissipation - of nanoscale systems. 
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